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ABSTRAGT 

We present polarization observations of two Galactic plane fields centered on Galactic coordinates 
(/,6) = (0°,0°) and (329°, 0°) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 
539 square degrees depending on frequency and field. These measurements were made with the 
QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of 
observations. The resulting maps represent the deepest large-area Galactic polarization observations 
published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 
2.8/iK deg, 2.3-6 times deeper than corresponding WMAP and Planck maps. The angular resolution 
is 27.'3 and 12.'8 EWHM at Q- and W-band, respectively. We find excellent agreement between the 
QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual 
instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from 
these in a manner consistent with reported systematic uncertainties for this channel. We combine 
QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining 
small angular scales from QUIET and large angular scales from WMAP. Erom these combined maps, 
we derive constraints on several important astrophysical quantities, including a robust detection of 
polarized synchrotron spectral index steepening of ^ 0.2 off the plane, as well as the Earaday rotation 
measure toward the Galactic center (RM= —4000 ± 200radm“^), all of which are consistent with 
previously published results. Both the raw QUIET and the co-added QUIET+ WMAP maps are made 
publicly available together with all necessary ancillary information. 

Subject headings: cosmic background radiation—cosmology: observations—Galaxy: general, center— 
polarization 
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1. INTRODUCTION 

The field of cosmic microwave background (CMB) cos¬ 
mology has undergone an important transition during 
the last two years. Until 2013, the primary limitation 
of most CMB experiments, whether targeting tempera¬ 
ture or polarization fluctuations, was instrumental noise. 
Contamination from astrophysical foregrounds and sys¬ 
tematic errors was generally small compared to the in¬ 
strumental noise level or intrinsic cosmic variance of the 
observations, and only minor corrections for either were 
required to produce robust cosmological results. Three 
examples among many are CBI (Readhead et al. 2004; 
Sievers et al. 2007), BOOMERanG (MacTavish et al. 
2006), and WMAP (Bennett et al. 2013). In this noise- 
dominated regime, the CMB field as a whole made rapid 
progress for more than two decades, with each new ex¬ 
periment improving cosmological parameter constraints 
in accordance with its respective noise level. Eventu¬ 
ally, this process led to the current highly successful 
ACDM ’concordance’ cosmological model, which today 
describes virtually all currently available cosmological 
observations with only six free parameters (Bennett et al. 
2013; Planck Collaboration XVI 2014; Planck Collabora¬ 
tion XIII 2015). 

This situation changed dramatically with the Planck 
release in 2013 March (Planck Collaboration I 2014), 
and later with the BICEP2 release the following year 
(BICEP2 Collaboration et al. 2014). The exquisite in¬ 
strumental sensitivity of Planck resulted in a CMB tem¬ 
perature likelihood that is, for the first time, limited by 
confusion from astrophysical foregrounds rather than in¬ 
strumental noise (Planck Collaboration XV 2014; Planck 
Collaboration XI 2015). Likewise, BICEP2 was the 
first CMB B-mode polarization experiment to become 
foreground-limited in polarization (BICEP2/Keck Array 
and Planck Collaborations 2015). 

To continue rapid progress towards a more refined 
cosmological model, in particular with respect to large- 
scale polarization, reionization, and inflation (e.g.. Pid¬ 
dle & Lyth 2000, and references therein), a thorough 
understanding of relevant astrophysical foregrounds is 
paramount. Great progress has already been made on 
this (e.g., Einkbeiner et al. 1999; de Oliveira-Gosta et al. 
2008; Bennett et al. 2013; Ichiki 2014), and in early 2015 
the Planck collaboration presented the most detailed 
full-sky model for the frequency range between 30 and 
353 GHz to date, including both polarized synchrotron 
and thermal-dust emission over the full sky (Planck Gol- 
laboration X 2015). According to this model, the fre¬ 
quency minimum for polarized foregrounds on degree an¬ 
gular scales occurs between 70 and 80 GHz, varying only 
weakly with multipole range, probably depending some¬ 
what on sky location. 

In order to improve on this foreground model, better 
measurements are required with respect to both depth 
and frequency coverage. In addition, control of instru¬ 
mental systematic errors is of course critical. As de¬ 
scribed in Planck Collaboration I (2015); Planck Col¬ 
laboration H (2015); Planck Collaboration VHI (2015); 
Planck Collaboration X (2015), there are still outstand¬ 
ing issues with the most recent Planck polarization obser¬ 
vations, both below and above the foreground frequency 
minimum at 70 GHz. Cross-checks and comparisons with 


external data sets, including WMAP^ can be helpful in 
identifying such issues. Other data sets anticipated in 
the very near future that should be useful in the effort 
to map out the foregrounds include S-PASS (2.3 GHz; 
Carretti et al. 2009), C-BASS (5 GHz; King et al. 2010, 
2014) and QUIJOTE (10-40 GHz; Rubiho-Martm et al. 
2012), all observing at low frequencies. 

In this paper, we present data that fit naturally into 
this larger astrophysical foreground program: measure¬ 
ments at 43 and 95 GHz of two fields in the Galactic plane 
taken by the QUIET instrument (QUIET Collaboration 
et al. 2012b) between 2008 October and 2010 Decem¬ 
ber. QUIET was a pathfinder experiment designed to 
improve limits on B-mode polarization and demonstrate 
the low level of systematic error achievable through the 
combination of careful monolithic-microwave-integrated- 
circuit (MMIC) receiver module design, instrument de¬ 
sign, and survey strategy. The instrument employed de¬ 
tector arrays comprising 19 Q-band (43 GHz) and 90 W- 
band (95 GHz) detector modules, observing from the At¬ 
acama Desert in Chile. The experiment reported the 
cleanest microwave polarization spectra with respect to 
instrumental systematic errors at the time. The sum of 
all instrumental systematic errors was constrained to cor¬ 
respond to a tensor-to-scalar ratio of r < 0.01 (QUIET 
Collaboration et al. 2012b). This result was only barely 
surpassed by the very recent and vastly more sensitive 
BICEP2 observations, which reported an equivalent limit 
on instrumental systematic errors of r < 0.006 (BICEP2 
Collaboration et al. 2015). Cosmological CMB E and B 
angular power spectra were reported in QUIET Collab¬ 
oration et al. (2011, 2012a), while constraints on polar¬ 
ized point sources were reported by QUIET Collabora¬ 
tion et al. (2014). 

The rest of this paper is organized as follows. In § 2 we 
review the QUIET data selection and processing pipeline 
as applied to the Galactic plane analysis, emphasizing 
those steps that are different compared to the original 
CMB-oriented analysis. We discuss the Q-band maps de¬ 
rived for the Galactic center field in § 3, while equivalent 
discussions and figures for the remaining observations are 
deferred to Appendix A. In § 4 we derive constraints on 
important astrophysical quantities such as the spectral 
index of synchrotron emission and the Earaday rotation 
measure toward the Galactic center, both of which are 
critical for performing robust astrophysical component 
separation. We summarize and conclude in § 5. All final 
data products (sky maps, mask, noise covariance matri¬ 
ces and beam profiles) are available on the LAMBDA 
website^^. Eollowing both WMAP and Planck^ we adopt 
the HEALPix (Gorski et al. 2005) convention for polar¬ 
ization, which differs from the lAU convention in the 
sign of the Stokes U parameter. All maps are provided 
in Galactic coordinates. 

2. OBSERVATIONS AND DATA PROCESSING 

The basic data selection and processing steps were de¬ 
scribed in detail in QUIET Collaboration et al. (2011, 
2012a). We briefly review the main steps in the following, 
describing in greater detail a few notable differences be¬ 
tween the previous and the updated analysis. The most 
important of these is co-addition with a second experi- 

http://lambda.gsfc.nasa.gov 
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TABLE 1 

Summary and comparison of field characteristics^ 



Q-band 

W-band 


Feature 

G-I 

G-2 

G-1 

G-2 

Field center, (/, 6) . 

(329°,0°) 

(0°,0°) 

(329°,0°) 

o 

O 

o 

sky area. 

483 deg^ 

301 deg^ 

573 deg^ 

539 deg^ 

Afpix (iV.ide = 512) . . 

36 831 

22 983 

43 668 

41090 

FWHM angular resolution 





QUIET. 

27.'3 


12.'8 


WMAP . 

30.'6 


13.'2 


Planck . 

27.'0 




Effective frequency, i/eff 





QUIET. 

43.1GHz 

94.5 GHz 


WMAP . 

40.5 GHz 

94.2 GHz 


Planck . 

44.1GHz 



Noise Q/U rms per 7' pixel 





QUIET. 

17 fiK 

2AiiK 

15 yK 

21 yK 

WMAP . 

58)uK 

64/rK 

96 yK 

108 yK 

Planck . 

52 yK 

55 yK 



Linear regression (y = ax A- h) slope, a 




(§ 3.2) 





a; =QUIET; y =WMAP, Q 

1.06 ±0.04 

1.05 ±0.04 



U 

I.IIT0.17 

1.00 ±0.03 



X =QUIET; y =Planck, Q 

1.33 ±0.30 

0.95 ±0.11 



U 

0.86 ±0.12 

1.00 ±0.04 



X = WMAP; y = Planck^ Q 

I.I9±0.27 

0.90 ±0.14 



U 

0.78 ±0.14 

1.00 ±0.05 



Noise-weighted mean and standard de¬ 
viation of deck split null map (§3.1) 





Stokes Q . 

0.22 ± 1.22/rK - 

0.32 ± 1.08 )UK 

0.09 ± 1.03 )UK -0.17 ± 1.01 )UK 

U . 

0.03±I.I6)uK 

0.10 ± 1.05 )uK 

-0.12 ± 1.04 yK 0.20 ± 1.01 yK 

QUIET systematic uncertainties 




(QUIET Collaboration et al. 2012b) 





Absolute responsivity . . . 

6% 


8% 


Absolute detector angle . 

I?7 


0?5 



a Note that Planck has not yet released W-band polarization maps, and the corresponding table entries are 
therefore empty. 


ment (typically WMAP), which is essential in order to 
robustly measure angular scales comparable to the size 
of the observed field. While CMB power spectrum or 
cosmological parameter estimation can be done without 
them, these modes are essential for deriving astrophysi- 
cal spectral parameters, such as the synchrotron spectral 
index or thermal dust temperature. A second difference 
is in the data-selection criteria, which are slightly less 
stringent in this Galactic analysis than in the original 
GMB analysis because the signal-to-noise ratio is higher 
for the Galactic fields. 

In the original QUIET analysis work, two pipelines 
were developed independently for cross-validation pur¬ 
poses, one implementing a pseudo-C^ framework (Hivon 
et al. 2002; QUIET Gollaboration et al. 2011), and the 
other implementing a brute-force maximum-likelihood 
framework (Tegmark 1997; Bond et al. 1998; QUIET Gol¬ 
laboration et al. 2011). A major advantage of the latter is 
that it provides both unbiased sky maps and correspond¬ 
ing dense pixel-pixel covariance matrices, which are use¬ 
ful for comparison and inverse-noise-variance co-addition 
with external data sets. The following analysis is based 
entirely on the maximum-likelihood pipeline. 


2.1. Observations and data selection 


The full, unfiltered data set consists of more than 
10,000 hours of observations taken from the Chajnantor 
plateau in Ghile between 2008 October and 2010 Decem¬ 
ber, covering two frequency bands (Q and W, with center 
frequencies 43 and 95 GHz, respectively) and six separate 
fields on the sky, not counting various calibration targets. 
Eour of these fields were selected as the primary low- 
foreground patches from which CMB constraints were de¬ 
rived (QUIET Collaboration et al. 2011, 2012a). To con¬ 
strain polarized foregrounds at the same frequencies, two 
high-foreground fields targeting the Galactic plane were 
also included in the observation schedule, and these mea¬ 
surements are the subject of this paper. These fields were 
observed when no primary CMB targets were available, 
effectively filling in right ascension “gaps” in the observ¬ 
ing schedule. Their positions on the sky are indicated in 
Eigure 1, together with corresponding patches of the po¬ 
larized synchrotron and thermal-dust maps recently pub¬ 
lished by Planck (Planck Collaboration X 2015). Adopt¬ 
ing the notation introduced in QUIET Collaboration 
et al. (2011), we will refer to the two Galactic fields as 
G-1 and G-2. The line-of-sight in G-1 cuts through the 
Centaurus arm as well as the tangent region of the Norma 
spiral arm (Garcfa et al. 2014), while G-2 contains the 
Galactic center region. Thus, these fields cover the two 
most populated areas of the Galactic disk in terms of 
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Fig. 1. — Positions of the QUIET fields superimposed on Planck 
foreground maps. The white and black outlines of the Galactic 
fields show the masks applied to the W-band (95 GHz) and Q- 
band (43 GHz) data, respectively. The central map is the Planck 
44 GHz temperature map, smoothed with a 1° FWHM Gaussian 
beam, with intensity range from -0.2 mK to 1 mK. The upper and 
lower panels compare the QUIET fields to the Planck polarized 
thermal dust and synchrotron foreground maps. Grid cell width is 
10 °. 

molecular gas (and therefore dust). 

Summary characteristics are provided for both fields in 
Table 1, including positions, sky area, effective noise lev¬ 
els and basic data-quality statistics. Regarding system¬ 
atic uncertainties, we include only the absolute respon- 
sivity and polarization angle in Table 1, as these are the 
most important ones for foreground analysis. We discuss 
the level of temperature-to-polarization leakage in our 
maps in § 3.1. For a complete discussion of systematic 
uncertainties relevant for B-mode analysis, we refer the 
interested reader to QUIET Collaboration et al. (2012b). 

The basic observing block of the QUIET scanning 
strategy was a so-called constant-elevation scan (CES), 
in which the sky drifts through the field-of-view while 
the telescope performs a simple periodic azimuth slew of 
fixed amplitude. For the four CMB fields and the G-1 
field, the scan width was 15°, while for G-2 it was 10°. 
Once the target field center drifted through the field-of- 
view by the same amount, the observing elevation was 
changed either to the opposite edge of the same field, 
thereby tracking the field on the sky in a set of discrete 
elevation steps, or to another field. The typical duration 
of a single CES was 30-60 minutes, depending on field 
size and elevation. In addition, the instrument was ro¬ 
tated about the optical axis (we refer to this as a deck 



-2000 -1000 0 1000 2000 -2000 -1000 0 1000 2000 

faK fiK 

Fig. 2.— Example of scan cut due to excessive striping. Both 
maps show the data taken by a single detector module during a 
single CES. The left-hand map shows a normal CES (CES id no. 
1808), while the right-hand map exhibits striping aligned with the 
scanning strategy (CES id no. 1826). All CES containing at least 
one such feature were cut from the analysis. Both example CES 
are taken from the W-band (95 GHz) observations of field G-1. 

rotation) in steps of 45°, typically once a week, for a to¬ 
tal of 8 angles. The combination of natural cross-linking 
from sky rotation and frequent deck rotations provided 
excellent modulation properties, suppressing many im¬ 
portant instrumental systematic effects (QUIET Collab¬ 
oration et al. 2012b). 

As described in QUIET Collaboration et al. (2011, 
2012a), the CMB-oriented QUIET data-reduction pro¬ 
cess was based on a fully-blind analysis philosophy, in 
which all data cuts, filters, and main processing steps 
were defined and finalized before inspecting any final sci¬ 
ence results, including power spectra and cosmological 
parameters. In this process, each CES-diode (i.e., the 
CES time-stream from an individual detector diode) was 
scanned for signs of contamination, and removed from 
the data set if any problematic issues were identified. 
The selection criteria assumed a low signal-to-noise ratio 
in any single CES, and relied strongly on accurate noise 
and modeling. This assumption, while valid for ob¬ 
servations of the low-foreground CMB sky, does not hold 
for the Galactic patches considered in this paper. Rather, 
the amplitude of the Galactic polarization signal is suffi¬ 
ciently high that the strongest signals are visible even in 
a single GES, in particular at 43 GHz. Under the default 
GMB-targeted selection criteria, a large fraction of the 
Galactic observations would be automatically excluded 
for this reason alone, and the default pipeline is there¬ 
fore not suitable for Galactic fields. On the other hand, 
the same high signal-to-noise ratio also implies that the 
fractional errors due to noise-modeling errors are much 
less important for Galactic than for CMB analysis. 

For these reasons, we adopt the following simplified 
data-selection criteria in this paper: First, we apply 
the default selection pipeline to eliminate obviously void 
scans, but exclude any tests that depend directly on fits 
of noise quality. The CES removed in this step include 
those affected by faulty hardware, and those for which 
the Moon was within the telescope’s sidelobes. Second, 
we manually remove scans with poor pointing (i.e., scans 
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TABLE 2 

Summary of data selection^ 


Q-band W-band 


Description 

G-1 

G-2 

G-1 

G-2 

Total CES count. 

295 

189 

568 

352 

Automatic cuts. 

1 

13 

27 

8 

Poor pointing. 

8 

0 

4 

4 

Short duration. 

3 

0 

14 

0 

Excessive striping. 

5 

2 

62 

27 

CES count after cuts . . . 

278 

174 

461 

313 

Observing time. 

. 303 h 

92 h 

595 h 

273 h 

efficiency . . . . 

94% 

91% 

83% 

90% 


“ Table lists the number of CES cut from the final data set by the 

cut criteria described in § 2.1. 

that did not hit the main target field) or short duration. 
Likewise we remove a small, discrete set of scans that on 
visual inspection was found to exhibit large-amplitude 
striping aligned with the scanning strategy. The cause 
of this latter effect is unknown. It is illustrated in Fig¬ 
ure 2 through two single-GES, single-detector-module 
maps from the W-band observations of field G-1. The 
left-hand panel shows a normal map, while the right- 
hand panel shows a contaminated map. The data con¬ 
taining the highest proportion of these GES are the W- 
band observations of field G-1, in which they make up 
10 % of the full data set. Although a substantial fraction 
of these could be saved through more aggressive filter¬ 
ing, considering the high signal-to-noise ratio of these 
Galactic sky maps we prefer to minimize systematic ef¬ 
fects rather than instrumental noise, and conservatively 
remove all GES that contain at least one striped single¬ 
detector map. 

Table 2 summarizes the data-selection statistics, both 
in terms of individual cuts and total observing efficiency. 
In total, 392 (868) observation hours are included in 
the final Q-band (W-band) maps, corresponding to an 
acceptance rate of about 90%, depending on field and 
frequency. For comparison, the acceptance rate ob¬ 
tained in the GMB-oriented QUIET analyses was 70- 
73% (QUIET Gollaboration et al. 2011, 2012a). 

2.2. Mapmaking 

Given a set of time-ordered data, we employ full 
maximum-likelihood mapmaking to estimate unbiased 
sky maps, m, by solving the normal equations (e.g., 
Tegmark 1997; Keskitalo et al. 2010) 

m=(p'^N“Tp) T'^N“Td. (1) 

Here symbols marked by denote pure time-domain ob¬ 
jects, while P and m denote (at least partially) map- 
domain objects. Specifically, P represents the pointing 
matrix, as defined by the default QUIET pointing model 
(QUIET Gollaboration et al. 2012a); N is the time- 
domain noise covariance matrix, given by a 1/f noise 
model; F represents a general time-domain filter; and d 
denotes the actual time-ordered data. The corresponding 
map-domain noise covariance matrix is 

N = (p'^N“Tp) ^ (p'^F'^N“Tp) (p'^N“Tp) ^. 

(2) 


Gonversion between measured voltages and thermody¬ 
namic temperature units, bandpass integration, and 
polarization-angle definitions are all implicitly described 
by the pointing matrix, P. For full details and conven¬ 
tions, see QUIET Gollaboration et al. (2011, 2012a,b). 

We use a HEALPix^^ Aside = 512 grid with 1' x 1' 
pixels for our maps, sufficient to support the 12.'8 FWHM 
angular resolution of the QUIET W-band channel. The 
QUIET Q-band channel has angular resolution of 27.'3 
FWHM, and could in principle be pixelized with 14' x 14' 
pixels; however, for consistency we pixelize both channels 
with the same grid. The WMAP polarization sky maps 
also use Aside = 512 grids. 

The total number of observed Q-band (W-band) pixels 
is 47,288 (52,391) and 45,162 (56,216) for G-1 and G-2, 
respectively. However, many of these pixels are observed 
only a few times, and therefore have low signal-to-noise 
ratio. In order to reduce the cost of subsequent matrix 
operations and data volumes, and acknowledging the fact 
that we later will co-add our maps with WMAP maps, 
we apply a mask to each full map, removing any pixels 
for which the effective QUIET noise rms is more than 
1.5 times the corresponding WMAP noise rms. The re¬ 
sulting masks are shown as black and white contours in 
Figure 1 for both Q- and W-band. 

While the QUIET GMB analysis used several comple¬ 
mentary time-domain filters, the requirements for fore¬ 
ground observations are somewhat different (see above). 
For these new maps, we have found that a mildly- 
apodized high-pass filter with a cut-off frequency of 
0.5 /scan is sufficient to produce maps with no obvious ar¬ 
tifacts, where /scan ~ 0.1 Hz is the scanning frequency of 
the telescope. The only exception is a specific large-scale 
feature in the G-1 field due to poor cross-linking. How¬ 
ever, as described below, rather than imposing a stronger 
time-domain filter in this case, we project out all har¬ 
monic modes with ^ < 10 in the pixel domain, to avoid 
excessive anisotropic filtering. No low-pass filters are ap¬ 
plied, in view of the fact that Galactic features tend to 
be strongly localized and full angular resolution is par¬ 
ticularly important. 

2.3. Co-addition with external data sets 

QUIET is for all practical purposes insensitive to phys¬ 
ical modes with wavelengths comparable to the size of the 
observed field. The QUIET field diameter of A ^ 20° 
thus suggests a loss of sensitivity for modes of ^ < 18^^. 
Although not vital for GMB power spectrum estimation, 
these modes are important for astrophysical foreground 
inference. We therefore co-add QUIET with an external 
large-scale experiment in order to produce optimal all¬ 
scale maps. Algorithmically, the co-addition is given by 
an inverse-variance-weighted sum of the form 



where the sums run over experiments, and represents 
the noise covariance matrix for the i’th experiment. The 

http://healpix.sourceforge.org 

In the QUIET CMB analyses, a lower limit of ^ = 25 was 
chosen for CMB power spectrum estimation from QUIET data 
(QUIET Collaboration et al. 2011, 2012a). 
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covariance matrix of the final map is 



( 4 ) 


If some set of N modes (which may be organized 
column-wise into an A^pix x N matrix V) happens to be 
affected by instrumental systematic errors in a given ex¬ 
periment, they can be projected out from the correspond¬ 
ing experiment covariance matrix before co-addition. 
This is most easily done by means of the Sherman- 
Morrison-Woodbury formula (e.g., Woodbury 1950), 

Nj Ni - NjV (V'^N“ W) V'^N*. (5) 

Effectively, this operation assigns infinite variance to all 
modes in V, ensuring that those modes do not contribute 
to the final map. In practice, we will use this operation 
to project out the largest-scale modes to which we can 
be certain that QUIET, due to its finite field size, has no 
sensitivity, by letting V consist of all spherical harmonics 
with i < 10. 

We use the notation described above to define an 
instrument-specific weight operator, F^, of the form 



which simply measures the relative weight carried by ex¬ 
periment i of each mode in the final map. Eor instance, 
Fquiq is the contribution from QUIET to the total map, 
mtot- Note that the sum over these operators is unity, 
= 1 , ensuring that the final map will be unbi¬ 
ased irrespective of instrument-specific filtering, as long 
as each individual map is inherently unbiased. 

Although QUIET, WMAP^ and Planck all nominally 
observe at Q-band, they do have slightly different band- 
passes and effective frequencies, as listed in Table 1. To 
account for these differences, we rescale the WMAP and 
Planck maps to the nominal QUIET frequency uq before 
co-addition, assuming a synchrotron-type power-law in¬ 
dex across the bands. Explicitly, the scaling factor for 
converting a map from frequency to is 


where 




5(0 


(e^ - 1)^ 

x‘^e^ 


hu 

^bTcmb 


( 7 ) 

(8) 


is the conversion factor between brightness and differen¬ 
tial thermodynamic temperature. Here, h and kB de¬ 
note the Planck and Boltzmann constants, and Tcmb = 
2.7255 K is the CMB monopole temperature. In order 
to avoid circularity in the analysis process, we adopt the 
synchrotron spectral index values reported by Euskeland 
et al. (2014) for these re-scaling factors, not those that 
will be derived from the QUIET maps themselves in § 4. 
Specifically, Euskeland et al. partitioned the whole sky 
(excluding bright compact objects and a region around 
the Galactic center of radius 1°) into 24 regions, and 


estimated the synchrotron spectral index for each re¬ 
gion from the WMAP K and Ka-band polarization sky 
maps. Eor G-1, we adopt the mean of their regions 23 
and 24 (see Eigure 1 in Euskeland et al. 2014), resulting 
in Pq.i = —2.93 ±0.01, while for G-2 we use the mean of 
regions 15 and 24, resulting in Pq _2 = —3.00 ±0.009. Eor 
WMAP^ these spectral indices translate into scaling fac¬ 
tors for G-1 and G-2 of 0.84 and 0.83, respectively. Eor 
Planck the corresponding factors are 1.067 and 1.069. If 
we instead were to adopt the spectral indices derived in 
§ 4 from QUIET (i.e., f3 = —3.12 ± 0.06), these num¬ 
bers would change by 1.2% and 0.7% for WMAP^ and 
by 0.4 % and 0.3 % for Planck. The impact of the precise 
value of the assumed spectral index is small compared 
to the intrinsic absolute responsivity uncertainty of 6 % 
in the QUIET observations (QUIET Gollaboration et al. 
2011). Eor W-band, the difference between the WMAP 
and QUIET frequencies is negligible, and we omit any 
re-scaling in this case. A Planck W-band polarization 
map is not yet available (Planck Gollaboration I 2015). 

2.4. Detailed analysis of Q-band G-1 field 

Before presenting the results from our analysis, it is 
useful to gain some intuition regarding the co-addition 
process described above. We therefore present the com¬ 
bination of the QUIET and WMAP Q-band G-1 maps 
in Eigure 3. The top row shows the raw maps from each 
experiment separately. The QUIET map is dominated 
by an essentially unconstrained mode with a gradient 
extending from the upper right to lower left. Qualita¬ 
tively similar features appear in all QUIET maximum- 
likelihood maps, but with an amplitude that varies 
strongly from field to field. In terms of how well the 
largest scales are constrained, the G-1 field is by far the 
worst of all six QUIET fields, whereas G-2 is one of the 
best. The fundamental difference between these two ex¬ 
treme cases lies in their degree of cross-linking coupled to 
the size and shape of the field. G-1 is neither a primary 
GMB field nor a particularly useful calibration source. 
As a result, it was observed through a relatively small 
set of scanning directions, from upper right to lower left 
edge in Eigure 3. Moreover, only half the amplitude 
was scanned within a single GES, resulting in two only 
partially overlapping and almost independent GES sets. 
Gonsequently, large-scale modes aligned with this direc¬ 
tion are poorly constrained. 

The G-2 field, on the other hand, covers the Galactic 
center, including the Galactic center arc, the strongest 
polarized object within the QUIET patches, and there¬ 
fore is a particularly useful calibration source, both for 
absolute responsivity and for pointing reconstruction 
(QUIET Gollaboration et al. 2012b). As a result, this 
field was observed from many different angles, leading to 
a more symmetric map. Additionally, the G-2 field was 
smaller to reduce the noise per sky area, and could be 
scanned from edge to edge within a single GES. Thus, 
even the large-scale modes are quite well constrained in 
G-2, and the raw G-2 map shows only weak evidence for 
spurious large-scale gradients. 

Substantial benefits can be obtained by co-adding data 
from QUIET with some large-scale experiment such as 
WMAP and/or Planck. The middle row of Eigure 3 
shows the inverse-variance-weighted component maps for 
QUIET and WMAP, Fqiiiq and Fw^w, as estimated 
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Fig. 3. — Example of map co-addition, applied to the Q-band 
QUIET and WMAP Stokes Q maps of the G-1 field. The top 
row shows the QUIET maximum-likelihood map, mg, and the 
WMAP 9-yr map, mw The middle row shows the contribution 
to the co-added map from each data set, computed by applying 
the corresponding weight operators, defined in § 2.3, to each map: 

Fq^niQ and F^^mw The co-added map, shown in the bottom 
left panel, is the sum of the two contributions. The framed panel 
(bottom right) shows an equivalent co-added map made using a 
version of the QUIET map where all modes of ^ < 10 have been 
discarded prior to co-addition; see § 2.4 for further details. 


from Equation 6. QUIET dominates the small-scale 
features in terms of signal-to-noise ratio, while WMAP 
dominates the large-scale modes. The previously domi¬ 
nating gradient in the raw QUIET map is effectively sup¬ 


pressed, and the weighted sum of the two contributions 
is shown in the bottom left panel. 

Nevertheless, a non-zero gradient is present even in the 
co-added map, albeit at a greatly reduced level. Neither 
WMAP nor Planck observe this mode; it is clearly due to 
an instrumental effect in QUIET, perhaps ground pick¬ 
up (QUIET Collaboration et al. 2012b). Irrespective of 
its origin, but confident it is an artifact in the QUIET ob¬ 
servations, we apply an additional real-space filter that 
projects out all harmonic modes with < 10 from the 
QUIET map, as described by Equation 5. We therefore 
adopt these few modes entirely from WMAP^ rather than 
estimating them in terms of a weighted mean of the two 
experiments. The result is shown in the lower right panel 
of Eigure 3, and this map appears astrophysically more 
reasonable than the unfiltered version in the lower left 
panel. We evaluate the need for this filtering by compar¬ 
ing the rms of the QUIET maps before and after applying 
it. In Q-band, filtering decreases the rms of the G-1 maps 
by more than 50%, whereas the corresponding value for 
field G-2 is a mere 2%. Similar results are found for 
the W-band maps. Hence we conclude that such filter¬ 
ing is prudent in the case of field G-1, but unneccessary 
for the far better constrained case of field G-2. All co¬ 
added maps for G-1 presented in the following have been 
derived using QUIET maps pre-filtered in this way. 

3. SKY MAPS 

We are now ready to present the Galactic field sky 
maps as observed by QUIET. In order to avoid excessive 
repetition, we will focus our discussion around the Q- 
band G-2 field, i.e., around the Galactic center region 
at 43 GHz. Corresponding plots and discussions for the 
remaining three fields (G-1 at Q-band, and both G-1 and 
G-2 at W-band) are given in Appendix A. 

3.1. QUIET-only maps 

The third row of Eigure 4 shows the inverse-variance- 
filtered QUIET G-2 maps, Fq mg. As described in the 
next section, we choose for now to include QUIET and 
WMAP^ but not Planck, in the weighted sum defined 
by Equations 3-6. Thus, the modes that are weighted 
down by the Fg operator are those for which WMAP has 
lower instrumental noise than QUIET, as measured by 
the respective noise-covariance matrices. This translates 
to the large-scale modes, as illustrated in § 2.4. 

The bottom row shows corresponding (half-difference) 
null-maps. These are derived by first dividing the full 
time-ordered QUIET data set according to the angle 
of deck rotation, with one set consisting of data taken 
with deck angles 0°, 90°, 180°, and 270°, and the other 
with deck angles 45°, 135°, 225°, and 315°. Independent 
maps are made from each subset, which are then sub¬ 
tracted. In the absence of systematic errors, any such 
null-map should contain instrumental noise only, and, 
as already mentioned, the QUIET analysis is fundamen¬ 
tally dependent on understanding null-maps. In the orig¬ 
inal GMB-oriented analysis, more than 20 different data 
splits were included. In this paper we focus on the deck- 
angle split alone, because it is the most stringent test 
for Galactic fields. Noise-weighted mean and standard 
deviation values of the deck-split null-maps of all four 
fields, given in Table 1, show that these maps are consis¬ 
tent with the expected Gaussian distribution. The only 





QUIET Collaboration 



Fig. 4.— Inverse-noise-variance-weighted Q-band (43 GHz) maps of the G-2 field (centered at Galactic coordinates (/,6) = (0°,0°)), for 
Planck, WMAP, and QUIET. Columns show, from left to right, Stokes Q, Stokes U, and polarization amplitude P = ■ The 

top row shows the Planck map mp, filtered to only contain the small-scale modes observable by QUIET, as determined by the QUIET 
weighting operator Fq defined in § 2.3. Rows 2 and 3 show the equivalent maps for WMAP and QUIET, respectively. The bottom row 
shows the half-difference null maps of the deck-split QUIET data; see § 3.1 for full details. 
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Fig. 5. — I-to-Q/U leakage in the QUIET Q-band (43 GHz) G-2 
field (centered on Galactic coordinates {l,b) = (0°,0°)). The top 
panel shows the WMAP9 Q-band temperature data. The bottom 
panel shows the half-difference map (niQ^ — mQ^)/2 of QUIET data, 
split by deck-angle (rotation angle of optical axis), in polarization 
amplitude P = \JQ‘^ + U^. Both maps have been smoothed to 
40' FWHM. Grid cell width is 5°. 


significant excesses in the Q-band G-2 null-map are two 
small-scale features, one toward the very Galactic center, 
the other toward a compact object at Galactic coordi¬ 
nates (/,6) = (353?17, 0?76) that is identified as PGCSl 
030 G353.17+00.76 in the Planck Gatalogue of Gompact 


Sources (Planck Gollaboration VII 2011). 

The most likely explanation for these excesses is so- 
called temperature-to-polarization (I-to-Q/U) leakage. 
Each QUIET MMIG module contains four detector 
diodes, two measuring Q, two measuring U, as defined 
by the local detector coordinate system (QUIET Gol¬ 
laboration et al. 2012b). Based on sky-dips (i.e., el¬ 
evation nods designed to monitor relative gain varia¬ 
tions) and lunar and Galactic observations, QUIET Gol¬ 
laboration et al. (2011) found that the instantaneous 
temperature-to-polarization leakage for the Q-band de¬ 
tectors was about 1 % in Q and 0.2 % in U. Modulation 
by both sky and deck rotations effectively suppresses this 
effect in final maps. The deck-angle null-test shown in 
Eigure 4 therefore provides a very strict upper limit on 
the net final leakage. 

To quantify this effect more accurately, we compare 
the null-map polarization amplitude with the WMAP 
Q-band temperature map in Eigure 5, both smoothed 
to an effective resolution of 40' EWHM to reduce noise. 
Gomparing the two maps visually, the qualitative corre¬ 
lation between the polarization excess and the temper¬ 
ature signal is obvious. Eurthermore, we find that the 
peak value of the polarization amplitude in the null map 
at the Galactic center is about 40 /iK, while the corre¬ 
sponding peak temperature amplitude is 60 mK. Thus, 
the net I-to-Q/U leakage is about 0.07% in the deck-split 
null map. In terms of total net polarization amplitude, 
this deck-split leakage corresponds to less than 4% of 
the full polarization signal of the Galactic center source. 
Again, after averaging over all possible polarization de¬ 
tector angles, these numbers will be significantly lower 
in the final maps. 

3.2. Comparison with Planck and WMAP 

The top two rows in Eigure 4 show the Planck and 
WMAP maps, scaled to the QUIET frequency as per 
Equation 7, and filtered with the QUIET weight oper¬ 
ator, i.e.. Equip and Equiw By removing the same 
large-scale basis functions from each map, all three can 
be directly compared without confusion from poorly con¬ 
strained large-scale modes. A quantitative comparison 
between the filtered Q-band QUIET, WMAP and Planck 
maps is given in Table 1 in the form of best-fit linear re¬ 
gression slopes (Petrolini 2014); corresponding W-band 
results are not provided, due to the very low signal-to- 
noise ratio of the WMAP W-band sky map and non¬ 
availability of the Planck 100 GHz map. 

Several interesting features may be seen in these maps. 
Eirst and foremost, it is reassuring that all three exper¬ 
iments observe the same broad structures, namely the 
positive Galactic plane and bright negative Galactic cen¬ 
ter in Stokes Q, and the negative ’wings’ in Stokes U. 
However, there are noticeable differences as well, the 
most important of which is the much lower noise of the 
QUIET maps. While only broad features may be iden¬ 
tified in the Planck and WMAP maps, even beam-sized 
features may be picked out by eye in the QUIET map. 

A second important but more subtle difference is the 
apparent amplitude of the Galactic plane in Stokes Q. 
Both QUIET and WMAP appear to be slightly brighter 
than Planck in the Galactic plane. This effect is vi¬ 
sually more striking in Eigure 6, where we show pair¬ 
wise difference maps between all three experiments, all 
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Fig. 6. — Pairwise differences of Q-band (43 GHz) maps of the 
G-2 (Galactic center) field, downgraded to HEALPix Ngide = 64 to 
suppress noise. All maps are weighted using the QUIET weight op¬ 
erator Fq, defined in § 2.3, retaining only the small-scale modes ob¬ 
served by QUIET in the differenced maps. Rows show, from top to 
bottom, QUIET-VFMAP, QUIET-P/anc/c, and WMAP-Planck. 
Golumns show Stokes Q and Stokes U. 

repixelized on a 55' x 55' grid (A^side = 64) in order to 
suppress instrumental noise. From top to bottom, the 
three rows show QUIET—1TM74P, QUIET—P/anc/c, and 
WMAP—Planck. While QUIET and WMAP are consis¬ 
tent, Planck clearly shows a deficit in the Galactic plane 
compared to the other two experiments. 

Another illustration of the same effect is provided in 
Figure 7, which shows the Stokes Q differences with re¬ 
spect to QUIET of Planck (blue curves) and WMAP (red 
curves) as a function of Galactic longitude, averaged over 
a |6| < 1?5 latitude band around the Galactic plane. This 
corresponds to the three center-most rows of pixels in 
Figure 6, although the evaluation was in fact performed 
on the full-resolution maps. For comparison, we have 
also plotted the corresponding mean of the WMAP Q- 
band temperature map as dotted lines, but scaled with 
a factor of 0.002. This signal would be produced by 
temperature-to-polarization leakage of 0.2%. The thin 
colored lines correspond to multiplying the QUIET map 
by ±6 %, which is the QUIET Icr uncertainty in absolute 
calibration (QUIET Gollaboration et al. 2011). Finally, 
we have marked the Galactic center region within 1° (i.e., 
the bright blue object in the Galactic center seen in Fig¬ 
ure 4) with a gray band. As noted in § 2, we assume a 
synchrotron spectral index of /3 g- 2 = —3.00 for this field. 



Galactic longitude [deg] 

Fig. 7. — Latitude-averaged difference between QUIET and 
WMAP (red) and between QUIET and Planck (blue) Q-band 
(43 GHz) maps of field G-2 (Galactic center), in Stokes Q, eval¬ 
uated over a latitude band around the Galactic plane of \b\ < 1?5. 
All maps are filtered with the QUIET weight operator Fq, de¬ 
fined in § 2.3, retaining only the small-scale modes that are well 
constrained by QUIET. The colored regions indicate the absolute 
QUIET calibration uncertainty of ±6 %. The dashed lines show 
the latitude-band-averaged WMAP Q-band temperature ampli¬ 
tude scaled by a factor of 0.002, providing a rough template of 
0.2 % temperature-to-polarization leakage. The gray region marks 
an area in longitude ±1° around the Galactic center within which 
all results are dominated by uncertainties in the foreground spec¬ 
tral index. 

as estimated by Fuskeland et al. (2014). However, the 
bright center object is not included in this estimate, as 
its properties appear qualitatively different from the sur¬ 
rounding environment. In addition, the amplitude of this 
object is very large, reaching 2600 /xK at 30' scales, and 
any difference maps like those considered here are there¬ 
fore very sensitive to spectral index uncertainties. For ex¬ 
ample, if the spectral index is —2.7 rather than —3.0, the 
effective difference between QUIET and WMAP would 
be —40 /xK, similar to what is seen in Figure 7. For now, 
we therefore exclude the central |/| < 1° region from our 
discussion, noting that further spectral index estimation 
work is required before conclusions can be drawn for this 
object. 

Focusing on the remaining longitude region with 
|/| > 1° in Figure 7, we see again the good agreement be¬ 
tween QUIET and WMAP. The thick red line fluctuates 
around zero with an overall mean and standard deviation 
of — I ± 3 /xK. In § 3.1 we derived an upper limit on the 
polarization amplitude uncertainty from temperature-to- 
polarization leakage of ^ 4% in the QUIET deck-angle 
null-map, which is further suppressed in the full map due 
to cross-linking. The uncertainty in the QUIET maps 
due to such leakage is therefore well within the level in¬ 
dicated by the ±6% uncertainty shown in the plot. In 
sum, we do not find any evidence for significant resid¬ 
ual temperature-to-polarization leakage in the full signal 
maps, either in QUIET or WMAP. 

For Planck^ we see a systematic positive excess, con¬ 
sistent with Figure 6. In principle, this excess could be 
due to several different effects. However, its magnitude 
and spatial pattern suggest temperature-to-polarization 
leakage, as discussed extensively in Planck Gollabora¬ 
tion H (2015). Compare the QUIET—P/ancA; residual to 
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Fig. 8. — Final co-added QUIET+l^MAP Q-band (43 GHz) maps of field G-2 (centered at Galactic coordinates (/, 6) = (0°,0°)). Using 
the weight operators defined in § 2.3, these are expressed as Fqiiiq + Fwmw. Left and right panels show Stokes Q and U, respectively. 
Grid cell side length is 5°. 


the upper dotted line, which indicates the mean WMAP 
Q-band temperature signal as a function of longitude, 
scaled by a factor of 0.002. Note in particular the corre¬ 
lated structures between longitudes I ^ —5° and —10°. 
The dotted line provides an approximate template of the 
temperature-to-leakage signal. However, this template 
is only approximate, since the detailed leakage pattern 
will additionally depend on the Planck scanning strat¬ 
egy and detector orientation, effects not accounted for 
here. Note, though, that such scanning strategy modula¬ 
tion can only reduce the correlation between the observed 
residual and the simple leakage model, never enhance it. 
Additionally, these features can not be due to intrinsic 
spectral index variations (unlike the Galactic center), be¬ 
cause WMAP agree very well with QUIET in this region 
despite having a longer relative frequency lever arm than 
Planck. 

Residuals at this level are consistent with the uncer¬ 
tainties for temperature-to-polarization leakage given in 
Planck Collaboration H (2015). Planck Collaboration 
X (2015) gives explicit leakage corrections based on de¬ 
tailed astrophysical foreground modeling; however, those 
models are necessarily associated with significant un¬ 
certainties, as they depend sensitively on both instru¬ 
ment and foreground models, and in particular on the 
assumed bandpass properties of the instrument. The 
Planck 44 GHz polarization map fails a few null-tests 
(Planck Collaboration I 2015), and is therefore not used 
in the Planck 2015 CMB likelihood (Planck Collabora¬ 
tion XI 2015). Thus, the new high-sensitivity QUIET 
maps presented in this paper represent a unique oppor¬ 
tunity to improve the Planck leakage model in future 
analyses, by virtue of providing a clean and direct refer¬ 


ence in the region of the sky with the highest signal-to- 
noise ratio. 

3.3. Co-added sky maps 

Given the qualitative differences between the Planck 
and QUIET maps reported above, we co-add the QUIET 
and WMAP maps in the current set of released maps, but 
not the Planck maps. Co-addition with Planck^ and any 
other experiment observing the same field, can always be 
performed later. Similarly, we co-add with the WMAP 
Wl, W2, and W3 differencing assembly maps at W-band, 
but not W4, since this particular differencing assembly is 
known to have significantly worse noise properties than 
the other three channels (Bennett et al. 2013). 

Eigure 8 shows the final co-added Q-band G-2 map in 
Stokes Q and U. Comparing this to the filtered maps 
shown in Eigure 4, the most noticeable differences ap¬ 
pear, as expected, near the edges of the field, where the 
QUIET signal-to-noise ratio deteriorates. In the full co¬ 
added map, the Stokes Q amplitude remains high along 
the Galactic plane to the very edge, where it tapers off 
in the filtered version. Corresponding maps for the other 
data sets (Q-band G-1 and W-band G-1 and G-2) are 
shown in Eigure 16 in the appendix. 

4. ASTROPHYSIGAL IMPLIGATIONS 

The sky maps presented in § 3 offer a fresh view of 
astrophysical foregrounds at microwave frequencies. In 
this section, we compare the co-added QUIET maps to 
WMAP and Planck^ and estimate both the polarized syn¬ 
chrotron spectral index and the Earaday rotation mea¬ 
sure in the Galactic center. 
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Fig. 9.— Comparison between the co-added QUIET+VFMTP Q-band (43 GHz, top left) and W-band (95 GHz, bottom left) maps with 
the WMAP K-band (23 GHz, top right; Bennett et al. 2013) map and the Planck thermal dust model (evaluated at 353 GHz, bottom right; 
Planck Collaboration X 2015), for field G-2 (centered on Galactic coordinates (/, 5) = (0°,0°)). All plots are generated using the Line 
Integral Convolution algorithm (Cabral & Leedom 1993). The color scale indicates polarization amplitude, P = \/Q‘^ + 1/^, while the flow 
stripes indicate magnetic-field orientation, i.e., rotated 90° with respect to the local polarization orientation. 


4.1. Visual inspection of sky maps 

In Figures 9 and 10 we compare our maps with the 
WMAP K-band map (Bennett et al. 2013) and the 
Planck 2015 map of thermal dust polarization (Planck 
Collaboration X 2015). All maps have been smoothed 
with a 30' FWHM Gaussian kernel to reduce instru¬ 
mental noise. The map of thermal dust is evaluated 


at 353 GHz, whereas the effective frequency of the K- 
band map is 22.4 GHz for a synchrotron-like frequency 
spectrum (Page et al. 2003). The color scale indicates 
the polarization amplitude, P, while the flow pattern 
traces field lines rotated by 90° with respect to the mea¬ 
sured EVPA (electric vector position angle), correspond¬ 
ing roughly to the magnetic field direction (this corre¬ 
spondence is exact under the assumption of optically thin 
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Fig. 10.— Same as Figure 9, but for G-1. 

synchrotron radiation). All plots are generated using an 
implementation of the Line Integral Convolution algo¬ 
rithm (Cabral & Leedom 1993) called Alice, provided 
in the HEALPix package. 

Eigure 9 shows the Galactic center field G-2. The 
magnetic field structure in the QUIET Q-band (top left 
panel) shows correlation with that in the WMAP K-band 
map (top right panel), showing not only the high qual¬ 
ity of the maps, but also indicating that the Q-band sky 
is dominated by synchrotron emission. According to the 
Planck 2015 astrophysical baseline model summarized in 



0 1000 2000 3000 


fiK 

Eigure 51 of Planck Collaboration X (2015), synchrotron 
emission should dominate over thermal dust emission by 
about an order of magnitude at Q-band, and our mea¬ 
surements are visually consistent with this picture. 

Three individual features are noteworthy in these low- 
frequency maps. Eirst, as is well-known, (e.g.. Page et al. 
2007) the orientation of the magnetic field at low Galac¬ 
tic latitudes is parallel to the Galactic plane. Second, the 
magnetic field lines to the north of the Galactic center 
form a ’U’ shape with an opening angle of ^ 45°. These 
field lines correspond to the origin of the polarized fila- 
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Fig. 11.— Planck 2015 baseline astrophysical foreground recon¬ 
struction in intensity for the G-2 field, as estimated with Commander 
(Planck Collaboration X 2015). From top to bottom and left to 
right, the panels show 1) CMB, 2) synchrotron, 3) spinning dust, 
4) free-free, 5) thermal dust, and 6) CO J=2—>-1. The region inside 
the white boundary in the CMB map has been replaced with a con¬ 
strained Caussian realization as part of the Commander algorithm. 

merits IX and XIV identified by Vidal et al. (2014a) in 
the WMAP sky maps. Third, the magnetic field lines 
in the Galactic center are rotated by an angle of almost 
90° with respect to the Galactic plane in Q-band, while 
at K-band this angle is about 70°. This is the expected 
signature of Faraday rotation, and we discuss it quanti¬ 
tatively in § 4.4. 

In the QUIET W-band map (bottom left panel of Fig¬ 
ure 9), we also see that the magnetic field is well or¬ 
dered and parallel to the Galactic plane at low latitudes. 
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Fig. 12.— Same as Figure 11, but for field C-1. 

Comparison with the Planck thermal dust map (bottom 
right) suggests that the QUIET map is dominated by 
dust rather than by synchrotron emission. Based on this 
qualitative analysis alone, we conclude that the effective 
polarized foreground minimum must lie between 43 and 
95 GHz, well separated from either side. Overall, the 
QUIET measurements are in good agreement with the 
Planck determination that the foreground minimum lies 
between 70-80 GHz (Planck Collaboration X 2015). 

Figure 10 shows the corresponding information for G- 
I, the field centered on (/, 5) = (329°, 0°). Although both 
maps have lower signal-to-noise ratio than those of G-2, 
the same qualitative conclusions hold. The Q-band map 




























QUIET measurements of the Galactic plane at 43 and 95 GHz 


15 






_Q 

I 

a 


a 




50 100 150 

WMAP K-band amplitude [/xK] 

Fig. 13.— Scatter plots between WMAP K-band (23 GHz) and QUIET Q-band (43 GHz) polarization amplitudes for fields G-1 (left 
column) and G-2 (right column), considering separately on-plane (|6| < 2?5; top row) and off-plane {\b\ > 2?5; bottom row) pixels. The 
dashed red lines indicate the best-fit power-law fit to each data combination. 


is clearly dominated by synchrotron emission and the W- 
band map is dominated by thermal dust emission. 

Gomparing the QUIET W-band and the Planck ther¬ 
mal dust maps, we note the presence of a bright lo¬ 
calized feature in the former at Galactic coordinates 
(/,6) = (326°,—2°) with no counterpart in the thermal 
dust map. We identify this as the supernova remnant 
(SNR) G326.3-1.8 (MSH 15-56), as discussed by Green 
(2009). It is a composite SNR that contains a shell 
with a relatively steep radio spectrum and an interior 
plerion component with a flat spectrum (Dickel et al. 
2000; Weiler & Sramek 1988). Erom the K, Ka, and Q 
bands, we estimate its spectral index as (3 ~ —2.7, consis¬ 
tent with a (relatively flat spectrum) synchrotron source. 
This value is steeper than the spectrum of the SNR shell 
(/d = —2.34) measured by Dickel et al. (2000) between 
0.408 and 14.7 GHz. This new measurement by QUIET 
might indicate a steepening of the polarized spectrum 
with frequency. 

Eigures 11 and 12 show the Planck baseline tempera¬ 
ture reconstruction in our fields, allowing for direct com¬ 
parison between our maps and individual astrophysical 


temperature components. Erom these we see that the 
bright G-1 source is indeed recognized as a synchrotron 
emitter in the Planck model, with little or no counterpart 
in any other component. 

4.2. Spectral index of polarized emission 

In this section we determine the effective power-law 
index between WMAP K-band and QUIET Q-band, 
which are heavily dominated by synchrotron emission. 
All maps are smoothed to a common resolution of 1° 
EWHM, slightly larger than the WMAP K-band beam 
of 53'. All spectral index estimates are derived from 
the polarization amplitude, P = The as¬ 

sociated bias from instrumental noise is corrected us¬ 
ing the asymptotic estimator (Montier et al. 2015; Vidal 
et al. 2014b), which is a generalization of the estimator 
first proposed by Wardle & Kronberg (1974) for the case 
where the uncertainties erg and au are non-identical. 

Many studies have reported a significant steepening in 
the synchrotron spectral index at high Galactic latitudes 
compared to the Galactic plane (e.g., Kogut et al. 2007; 
Euskeland et al. 2014). As a first test, we therefore parti- 
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TABLE 3 

Polarized synchrotron spectral index between K and Q 

BAND^ 


Field 

QUIET 

WMAP 

G-1 full. 

-2.91 ±0.01 

-3.02 ±0.03 

|fe| < 2.5° . . . 

-2.89 ±0.04 

-3.02 ±0.06 

|6| > 2.5° . . . 

-3.14 ±0.03 

-3.17 ±0.08 

G-2 full. 

-3.01 ±0.01 

-3.02 ±0.02 

|6| < 2.5° . . . 

-2.89 ±0.02 

-2.93 ±0.03 

|6| > 2.5° . . . 

-3.07 ±0.04 

-3.16 ±0.07 


“ Spectral indices are computed between WMAP K-band and ei¬ 
ther QUIET (second column) or WMAP (third column) Q-band. 

tion our fields into high (|6| > 2.5°) and low (|6| < 2.5°) 
latitudes, and determine the spectral indices separately 
for each region. This calculation is summarized in the 
scatter plots of Figure 13 (based on polarization ampli¬ 
tudes) between WMAP K-band and QUIET Q-band for 
both fields (left and right columns) and low and high 
latitudes (top and bottom rows). The red dashed lines 
indicate the best-fit lines, corresponding to the spectral 
index indicated in each panel. We find results consistent 
with earlier measurements. In G-2 we derive a Galactic 
plane spectral index of [3 = —2.89 ± 0.02, which steep¬ 
ens to /3 = —3.07 ± 0.04 at high latitudes^^. For G-1, 
the corresponding numbers are (3 = —2.89 ± 0.04 and 
(3 = —3.14 ± 0.03. In both cases, the spectral index is 
about 0.2 steeper off the plane than in the plane. 

Including all latitudes in the spectral index evaluation, 
we find Pg-i = -2.91T0.01 and /3 g- 2 = -3.01T0.01, as 
listed in Table 3. These are fully consistent with the val¬ 
ues adopted from Fuskeland et al. (2014) (/dc-i = —2.93 
and /3 g- 2 = —3.00; see § 2.3), derived from WMAP K- 
and Ka-band, that were used to re-scale the WMAP and 
Planck maps to the effective QUIET Q-band frequency. 

4.3. Electric vector position angles (EVPA) 

As noted in § 4.1 and seen in Figures 9 and 10, the 
observed orientation of the polarization vector (EVPA) 
is roughly perpendicular to the Galactic plane in the mil¬ 
limeter wavelength range. Given the assumption of opti¬ 
cally thin synchrotron radiation, the Galactic magnetic 
field is therefore roughly parallel to the plane. To quan¬ 
tify these alignments, we measure the EVPA, 

X = iarctan(f7/Q)29 (9) 

defined as zero for polarization vectors aligned with 
Galactic meridians and increasing counter-clockwise in 
Galactic coordinates. Spatial variations in y can be used 
to constrain models of the Galactic magnetic field (e.g., 
Jaffe et al. 2013). 

In Table 4 we tabulate the mean EVPA as a function of 
Galactic latitude for both G-1 and G-2, averaging over 
latitude slabs of 2° width and masking bright objects. 
Averaging over the entire fields, the Q-band EVPA is 
-1?9 ± 0?4 (stat) ±1?7 (syst) for G-1 and -10?8 ± 0?6 
(stat) ±1?7 (syst) for G-2. Similar results were reported 

All error bars include the uncertainty in the QUIET absolute 
responsivity. 

29 yye remind the reader that QUIET uses the COSMO conven¬ 
tion for the Stokes parameters, where Stokes U is replaced by — U 
to follow the lAU convention. 


TABLE 4 

EVPA AS A FUNCTION OF GALACTIC LATITUDE^’^ 


h X in G-1 X in G-2 


Min Max Q-band W-band Q-band W-band 


+3? 

0 

±5? 

0 . . . 

30?0 

± 

0?3 

11?7 

± 

0?2 

3?5 

± 

0?3 

-20? 

7 

± 

0?2 

±1? 

0 

+3? 

0 . . . 

14?9 

± 

0?2 

7? 5 

± 

0?2 

0?3 

± 

0?1 

-17? 

2 

± 

0?2 

-1? 

0 

±1? 

0 . . . 

-0?8 

± 

0?1 

4?6 

± 

0?1 

7?0 

± 

0?1 

-4? 

7 

± 

0?1 

-3? 

0 

-1? 

0 . . . 

-5?3 

± 

0?2 

5?0 

± 

0?2 

10?1 

± 

0?2 

-12? 

9 

± 

0?2 

-5? 

0 

-3? 

0 . . . 

-9?9 

± 

0?4 

3?2 

± 

0?2 

-4?7 

± 

0?3 

-17? 

4 

± 

0?3 


X is measured counterclockwise relative to the Galactic north 
direction. 


^ Errors are statistical only. See Table 1 for systematic errors. 



ElC. 14.— Polarization intensity maps of the Galactic centre 
region seen by QUIET at 43 GHz (left) and 95 GHz {right), at the 
original angular resolution of 27^3 and 12'8 EWHM respectively. 
The vectors are rotated 90° with respect to the measured EVPA, 
indicating the orientation of the magnetic field. Both maps are 
displayed in Galactic coordinates with a grid size of 1°. 

by Bierman et al. (2011) for 100, 150, and 220 GHz de¬ 
rived from the BIGEP observations. 


4.4. Earaday rotation at the Galactic center 

The Galactic center patch G-2 includes the brightest 
polarized source observed by QUIET. Figure 14 shows 
the QUIET maps of the central region at 43 and 95 GHz. 
The magnetic field of the central source is well aligned 
and perpendicular to the Galactic plane. In the Galactic 
plane itself, the magnetic field is parallel to the plane 
in both bands. At Q-band, a ring with a minimum of 
polarized emission is seen around the central source. This 
is due to the cancellation of emission from the source and 
the Galactic plane with different polarization angles. 

The polarization of the Galactic center has been stud¬ 
ied in detail at higher angular resolution. Haynes et al. 
(1992) mapped the region at 8.5 GHz with a beam size 
^2f8 using the Parkes telescope. Tsuboi et al. (1995) also 
observed this region at 42.5 GHz (similar to QUIET-Q 
band) with 39" angular resolution. From these results, it 
is clear that the dominant source of polarized radio emis¬ 
sion originates from the Galactic center ‘arc’. This arc 
is a very narrow non-thermal filament, ~ 15' in length, 
perpendicular to the Galactic plane, located within a dis¬ 
tance of about 20' from the Sgr A* radio source at the 
Galactic center (see e.g. Yusef-Zadeh et al. 1984; Pedlar 
et al. 1989). 

We measure the EVPA of the Galactic center arc 
as a function of frequency from the QUIET, WMAP^ 
and Planck maps, all smoothed to a common resolution 
of 1° EWHM, averaging over a 20' diameter disk cen¬ 
tered on the peak of the polarized emission at (/, 6) = 
(0°15, —0°1). The resulting values are listed in Table 5, 
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TABLE 5 

Galactic center source polarization parameters^ 


Band 

Stokes Q 
[mK] 

Stokes U 
[mK] 

Pol angle 
[°] 

QUIET Q ... 

. -1.39 ±0.04 

0.01 ±0.02 

90.1 ±0.8 

W . . . 

. -0.89 ± 0.06 

0.67 ±0.03 

108.5 ± 1.3 

WMAP K ... 

. -0.44 ± 0.06 

-0.72 ±0.03 

60.5 ± 1.4 

Ka . . 

. -1.06 ±0.05 

-0.38 ±0.03 

80.1 ± 1.2 

Q ... 

. -1.39 ±0.04 

-0.17 ±0.03 

86.6 ±0.9 

V ... 

. -1.62 ±0.03 

0.26 ± 0.02 

94.6 ± 0.5 

W . . . 

. -0.84 ±0.06 

0.56 ± 0.04 

106.9 ± 1.6 


“ Values correspond to the mean evaluated inside a 20' diameter 
aperture centered at Galactic coordinates {l,b) = (0?15,—0?1). 
Uncertainties in Q and U correspond to the statistical fluctua¬ 
tions in the map, as measured in a larger aperture around the 
source. See Table 1 for systematic errors. 
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Fig. 15. — EVPA as a function of wavelength squared for the 
Galactic center arc for both WMAP (red points), and QUIET (blue 
points). The dashed line is a linear fit, as expected for a pure 
foreground Faraday screen resulting in an EVPA x oc RM A^. We 
exclude the W-band points from the fit; see § 4.4 for explanation. 
The best-fit rotation measure is RM = —4000 ± 200 rad m“^. 


including individual Stokes parameters and EVPAs. The 
latter are also plotted in Eigure 15 as a function of 
wavelength squared, A^. Due to Earaday rotation (e.g., 
Burn 1966), the observed EVPA is expected to follow 
X = RMA^, where RM is the rotation measure. 

Neglecting the W-band data points, we find a tight 
relation with the expected form, with a best-fit rota¬ 
tion measure of RM = —4000 ±200 rad m“^, correspond¬ 
ing to a total rotation angle of ~ 35° between 23 and 
60 GHz. Tsuboi et al. (1995) reported a value of RM 
= —3120 ± 188 rad m“^ between 10 GHz and 42.5 GHz 
within a 2' beam centered at the Galactic center arc. The 
difference between our value and the one from Tsuboi 
et al. (1995) could be due to Earaday depolarization at 
low frequencies, which biases the RM measurements to¬ 
wards lower values (Law et al. 2008). Large absolute 
values of RM similar to the one we reported here have 
also been measured close to the Galactic center by Roy 
et al. (2005). 

The EVPAs for the Galactic center arc found in the 
QUIET sky maps are in good agreement with those re¬ 
ported by WMAP^ both at Q- and W-band. However, 


the W-band angles differ significantly from the expected 
Earaday rotation prediction by about 15°. The reason for 
this is the contribution of dust polarization, which dom¬ 
inates W-band, while synchrotron only dominates below 
- 60 GHz. 


5. SUMMARY 

We have presented polarization measurements of the 
Galactic plane at 43 and 95 GHz, as observed with 
QUIET between 2008 October and 2010 December. The 
resulting Galactic maps are the deepest published to date 
at their respective frequencies, by a factor of 2-4 at Q- 
band compared to Planck and WMAP^ and by a factor 
of 5-6 at W-band compared to WMAP. We find no sig¬ 
nificant evidence of residual instrumental systematic er¬ 
rors in these high-signal-to-noise-ratio maps. We derive a 
conservative upper limit on temperature-to-polarization 
leakage of < 0.07 % in the Q-band, translating to a < 4 % 
uncertainty in polarization amplitude at the Galactic 
center. Eor comparison, the uncertainty in absolute re- 
sponsivity is 6 % for Q-band and 8 % for W-band. 

Our maps agree very well with corresponding WMAP 
polarization observations in both Q- and W-band over 
the entire fields. Accounting for the different effective 
frequencies and uncertainties in the synchrotron spec¬ 
tral index, we find no compelling evidence for signifi¬ 
cant systematic differences between the two. At 44 GHz, 
both QUIET and WMAP observe a stronger polariza¬ 
tion signal than Planck along the Galactic plane. One 
potential explanation for this difference is temperature- 
to-polarization leakage of ^ 0.2 % in the Planck data set, 
which possibly might be related to the null-test failures 
for this particular channel already reported by Planck 
Gollaboration H (2015). Eurther work is needed to un¬ 
derstand these discrepancies in detail. 

Exploiting the agreement between QUIET and 
WMAP^ we have computed inverse-noise-variance- 
weighted averages between the two experiments, com¬ 
bining small-scale information from QUIET with large- 
scale information from WMAP. The resulting maps are 
publicly available on LAMBDA, and should prove useful 
both for experimental consistency checks, as exemplified 
in this paper, and for understanding the physical proper¬ 
ties of polarized foreground emission at microwave wave¬ 
lengths. In the current paper, we have presented a few 
examples of such analyses, evaluating the spectral index 
of synchrotron emission, the mean EVPA near the Galac¬ 
tic plane, and the Earaday rotation measure of the Galac¬ 
tic center source. A key result from this work is robust 
evidence for true spatial variations in the synchrotron 
spectral index of diffuse polarized emission along and off 
the Galactic plane. 

Bruce Winstein, who led the QUIET project, died in 
2011, soon after observations concluded. The project’s 
success owes a great debt to his intellectual and scientific 
leadership. 

Support for the QUIET instrument and operation 
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PHY-0355328, AST-0448909, AST-1010016, and PHY- 
0551142; KAKENHI 20244041, 20740158, and 21111002; 
PRODEX G90284; a KIPAC Enterprise grant; and by 
the Strategic Alliance for the Implementation of New 
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APPENDIX 


SUPPORTING FIGURES 

In § 3 we presented both the raw QUIET and the co-added QUIET+IUMHP Q-band sky maps as derived for 
the Galactic center field, G-2. We also showed an internal consistency test for this field between QUIET, WMAP, 
and Planck^ in the form of difference maps. In this appendix we show corresponding plots for the remaining three 
data combinations, namely the Q-band G-1 and W-band G-1 and G-2 fields. All main conclusions remain unchanged 
compared to the original discussion, and the following plots are reproduced for completeness and reference purposes 
alone. Null-map statistics for all fields and and best-fit linear regression slopes for the Q-band data are listed in 
Table 1. 

Eigure 16 shows the final co-added QlJIETp WMAP Stokes Q and U parameter maps for Q-band G-1 (top row), 
W-band G-1 (middle row), and W-band G-2 (bottom row), corresponding to Eigure 8 for the Q-band G-2 field in the 
main text. Gomputing the polarization amplitudes and EVPAs from these leads to the maps shown in Eigures 9 and 
10 . 

Eigure 17 shows the (filtered) QUIET contributions to the co-added sky maps and the deck-split half-difference maps 
for each field combination, corresponding to Eigure 4 in the main text. No significant residuals are seen in any of 
these difference maps. The Q-band G-2 case discussed in the main text therefore represents a conservative worst-case 
scenario with respect to temperature-to-polarization leakage. 

Eigure 18 shows the difference maps between QUIET, WMAP^ and Planck^ all downgraded to A^side = 64 as in 
§3.2. Differences with respect to Planck are only evaluated for the Q-band G-1 field, since Planck does not provide a 
polarized W-band map at this time. As for the G-2 field, we note a significant positive residual with respect to Planck 
in G-1, while no significant residuals are seen with respect to WMAP in either case. Einally, it could be noted that 
the W-band difference maps have a higher noise level. This is due to the high noise rms in the WMAP W-band maps, 
as listed in Table 1. 
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Fig. 16.— Final co-added QUIET+VFMTP maps. The top row shows the Q-band (43 GHz) result for field G-1 (centered on Galactic 
coordinates (I, b) = (329°, 0°)). Using the weight operators defined in § 2.3, these maps are expressed as FqmQ +Fwmw- The middle and 
bottom rows show the equivalent maps for W-band (95 GHz), both fields. The left and right columns show Stokes Q and U, respectively. 
The grid cell width is 10°. 










QUIET measurements of the Galactic plane at 43 and 95 GHz 


21 


a 

E 

a 


Q-band, 

G-1 


-oQ 

E 


i-HlCsl 













. 







W-band, 


E 

I 


-oa 

1—l|CNI 


N I ^ 

-100 0 100 




fiK fiK fiK 

Fig. 17.— Inverse-noise-variance-weighted QUIET maps and deck-split half-difference maps. Columns show, from left to right, Stokes 
Q, Stokes U, and polarization amplitude P. The first row shows the QUIET map mg, for the G-1 field (centered on Galactic coordinates 
(/, b) = (329°, 0°)) in Q-band (43 GHz), filtered to only contain the small-scale modes observable by QUIET, as determined by the QUIET 
weighting operator Fq defined in § 2.3. The second row shows the corresponding deck-split half-difference map, |(mg^ — mg^). The 
following rows show the equivalent pairs of maps for both fields in W-band (95 GHz). 
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Fig. 18. — Pairwise differences of QUIET and WMAP or Planck maps, downgraded to HEALPix Ngidf. = 64 to suppress noise. All 
maps are weighted using the QUIET weight operator Fq, defined in § 2.3, retaining only the small-scale modes observed by QUIET in 
the differenced maps. The top row shows QUIET—VFMAP for field G-l (centered on Galactic coordinates {l,h) = (329°, 0°)) in Q-band 
(43 GHz). The second row shows the corresponding difference of QUIET—P/anc/c, while the third and fourth rows show QUIET—VFMAP 
for W-band (95 GHz), both fields. Golumns show, from left to right, Stokes Q and Stokes U. 






